Effects of Defects on Friction for a Xe Film Sliding on Ag(lll) 
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The effects of a step defect and a random array of point defects (such as vacancies or substitu- 
tional impurities) on the force of friction acting on a xenon monolayer film as it slides on a silver 
(111) substrate are studied by molecular dynamic simulations and compared with the results of 
lowest order perturbation theory in the substrate corrugation potential. For the case of a step, the 
magnitude and velocity dependence of the friction force are strongly dependent on the direction of 
sliding respect to the step and the corrugation strength. When the applied force F is perpendicular 
to the step, the film is pinned for F less than a critical force F c . Motion of the film along the 
step, however, is not pinned. Fluctuations in the sliding velocity in time provide evidence of both 
stick-slip motion and thermally activated creep. Simulations done with a substrate containing a 5 
percent concentration of random point defects for various directions of the applied force show that 
the film is pinned for the force below a critical value. The critical force, however, is still much lower 
than the effective inertial force exerted on the film by the oscillations of the substrate in experiments 
done with a quartz crystal microbalance (QCM). Lowest order perturbation theory in the substrate 
potential is shown to give results consistent with the simulations, and it is used to give a physical 
picture of what could be expected for real surfaces which contain many defects. 



I. INTRODUCTION 



In previous simulations of a xenon film sliding on a 
silver substrate, using a periodic (i.e., defect-free) sub- 
strate, a viscous force of friction (i.e. one proportional 
to the sliding velocity) was found ||,|, in agreement 
with the experimental results of Krim et. al. (^^]. In 
contrast, perturbation theory calculations give a veloc- 
ity independent contribution to the friction (i.e., "dry 
friction") when there are point defects (i.e., a point de- 
fect denotes a defect which is centered around a point 
in the lattice), such as vacancies or substitutional impu- 
rities [^-U in the substrate. Since, real surfaces, even 
very smooth ones, always contain defects, the viscous 
friction found in these experiments (3|J] is a surprising 
result. Perturbation theory results for straight line de- 
fects, i.e., defects that extend along a line such as steps or 
facet boundaries Q], however, are consistent with viscous 
friction. In this article, we perform molecular dynamics 
simulations for a film of Xe atoms sliding on a Ag(lll) 
substrate containing a step and also for a substrate with 
a random array of point defects. The point defects are 
found to pin the film for an applied force below a criti- 
cal value, whereas a step generally does not pin the film. 
The pinning force due to a 5 percent concentration of 
point defects comparable in strength to the corrugation 
potential, however, is still much lower than the effective 
inertial force per film atom in the experiment done by 
krim. |5|4|] this force can be expressed as muj 2 A, where 
m is the atomic mass and to and A are the frequency and 
amplitude, respectively, of the quartz crystal microbal- 
ance used in the experiment. Lowest order perturbation 
theory in the substrate corrugation potential for the sys- 
tem treated in the simulations is found to be consistent 



with the simulations. The perturbation theory results are 
then used to try to give a physical picture of what one 
would expect for real substrate surfaces, which contain 
many defects. 



II. SIMULATIONS 

A. The model used in the simulations 

The model Hamiltonian used in Ref. Q for N film 
atoms of mass m at positions r& (fc = 1, .., N) is given by 
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(2.1) 



where p^. is the momentum of the atom fc, and the total 
potential U{t\, • • • , rjv) is given by 

N N 

U(r 1 ,---,r N ) = '£u s (r k )+ £ VQ rj - r fe |). (2.2) 

k=l j<k=l 

Here, U s (rk) is a single particle potential describing the 
interaction between the fc-th film atom and the substrate, 
and V(\rj — r^l) is the pair potential interaction between 
the j-th and fc-th atoms in the film. 

The interaction between two Xe atoms is given by a 
Lennard-Jones potential 



V(r) = 4e 



(?)"-(?)' 



(2.3) 
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where e — 19.83 meV, and a = 4.055 A. The interaction 
between a Xe atom and the substrate can be described by 
a substrate potential without internal degrees of freedom 
given by § 



initially in the gas phase. The atoms condense in 250 io 
or less, forming a triangular lattice incommensurate with 
the substrate fcc(lll) surface. 



U s (r h z) = U (z) + U 1 (z)}2cos(G-r ll ), (2.4) 

{G} 

where rii = (x, y) are the coordinates of the Xe atom 
parallel to the substrate, and {G} is the set of the six 
shortest reciprocal latti ce vectors of the substrate. The 
first term in Eq. ( |2.4|) describes the mean interaction 
of the atoms with the substrate, and the second term 
describes the periodic corrugation potential. 

Expressions for Uq(z) and U\(z) were derived by Steele 
H assuming that the substrate potential U s (r) is a sum 
of Lennard-Jones potentials between one film atom and 
all of the atoms in the substrate. However, a potential 
like U s (r), which is a sum of Lennard-Jones potentials is 
not a correct description of the interaction of a metallic 
surface with a noble gas atom. The corrugation potential 
is reduced (from the value found by summing Lennard- 
Jones potentials) due to electronic screening. For this 
reason we employ a weaker corrugation potential, as did 
Cieplak et al. in Ref. |Q. The corrugation potential we 
use is 
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where a = iirexe/AgA* 6 



/V3, z* = z/a, a = 2.892 A 
is the lattice constant of the substrate, A* — (?xe/Ag / a i 
cji = 4n/^/3. We calculate the Lennard-Jones parameters 
axe/Ag and £xe/Ag by fitting (i) the position of the min- 
imum of Uq(z) to the distance between a Xe atom in the 
first layer and the ion cores of the substrate (zq ) , and (ii) 
the attractive well depth to the binding energy of one Xe 
atom to the Ag substrate (Uo(zo) — —211 meV, from fl). 
We find uxe/Ag = 4.463 A and e Xe /A g — 13.88 meV. 

The corrugation potential U\{z) in Eq. ( |2.5| ) falls off 
exponentially at large z. The above parameters give 
U\(zq) J2{g] cos (G ■ = 2.025 meV (for the maximum 
value of this sum) in contrast to the corresponding value 
of 10.13 meV for Steele's potential at z M. Our corru- 
gation gives good agreement with the experimental value 
of the slip time. 

Our simulations are carried out at an equilibrium tem- 
perature of T = 77 'A °K, and the particles move in a 
three dimensional box of size 20 ax lOa-s/3 x 10 a. The 
time scale for vibrations of the adsorbed film atoms is 
to = y/(ma 2 /e) = 3.345 ps , with m = 2.16 10~ 22 g. 
Periodic boundary conditions in the x and y directions 
are employed along with a hard wall boundary condition 
in the z direction at the top of the box. 

We change the coverage by changing the number of 
Xe atoms N. We use 60 < N < 370. All atoms are 



For the simulatio ns d one with a step present, the po- 
tential in equation (2.4) has the z-coordinate replaced by 
z—g{x). The use of this function guarantees that we have 
the same corrugation as we used for a substrate free of de- 
fects in Ref. @. Rereg(x) = 0.58a[f(x 1 -x)-f(x-x 2 )] 
where f(x) is the Fermi function, l/{e~ x /^ + 1]. (See 
Fig. [[].) We choose w to be equal to 1.1 a, AS THE 
width for the step edge, where a is the distance param- 
eter for the Lennard-Jones potential between a film and 
substrate atom. 




FIG. 1. A plot of g(x) = 0.58(7 [/(si - 
f(x) is the Fermi function. 



■x) — f(x— X2)], where 



We take x\ and X2, the locations of the beginning and 
end of the step, equal to L x /3 and 2L x /3, where L x is the 
length of the box (along x) in which the simulations are 
done. We thus assume a straight step in the substrate, 
(see Fig. |^ (a and b)), which runs along the y-axis. The 
height of the step is about 0.8 of an atomic distance. The 
use of a function g(x), which varies smoothly with x, is 
a reasonable choice because the nonzero radius of a sur- 
face atom makes the potential that acts on the atom vary 
smoothly as the atom moves over the surface. 



Most of our simulations were done with a coverage cor- 
responding to the uncompressed monolayer (163 parti- 
cles). Periodic boundary conditions in the x and y di- 
rections are employed along with a hard wall boundary 
condition in the z direction at the top of the box. Because 
of our use of periodic boundary conditions, we are tech- 
nically simulating a periodic array of defects. We feel, 
however, that the cell length used in our simulations of 
about 13 xenon atom spacings is sufficiently long for the 
coherence from one cell to the next to be unimportant 
for most values of the sliding velocity. This is justified 
using perturbation theory for reasonable values of 
the phonon damping constant in appendix A. 



2 



J" 




L,(o) 




FIG. 2. a. Upper view of the position of the particles after 
200,000 iterations of the program. The step is located along 
the y axis between 5 and 10 a. b. Side view of the same sets 
of positions. 

In the present simulations, we use the same method 
as Cieplak, et. al. jlj], in which a constant external force 
F is applied to each atom in the plane of sliding, and 
the resulting steady-state velocity is calculated (see Fig. 
|j) . Throughout the simulation a thermostat that rescales 
the three velocity components is used to maintain con- 
stant temperature. The rescaling is done in the center 
of mass reference frame so as not to change the cen- 
ter of mass velocity and thus introduce an unphysical 
force of friction due to the thermostat. It renormalizes 
the atomic velocities every time step, so that the total 
kinetic energy per atom in the center of mass frame is 
maintained at (3/2)fc^T, where fcs is Boltzmann's con- 
stant and T is the desired absolute temperature of the 
system (T = 77.4 "if in our case). This method allows 
us to determine conveniently the velocity dependence of 
the friction force. The velocity of the film is affected 
by the inhomogeneities of the medium, in this case the 
corrugation and the defects. In the absence of defects, 
the film can move for an arbitrarily weak applied force, 
but when there is a step present, the film can be pinned 
for F less than a critical value F c , if the force is applied 
perpendicular to the step. In other words, there are the 
following two phases, a pinned phase when the external 
force F is F < F c , and a moving phase. The depin- 
ning transition takes place at a critical threshold force 
F c . In the vicinity of the depinning transition the aver- 
age velocity has the form V ~ (F — F c )@ where (3 is the 
velocity exponent. 

The potential for a substrate containing point defects 
is of the form: 



V i (f-f f d ) = Voe 



(2.6) 



where Vb and I are the strength and range of the poten- 
tial respectively and is the location of the defect. The 
defect positions are chosen to be potential minima of the 
corrugation potential chosen at random with probability 
c, where c is the defect concentration. (The simulations 
done for point defects were done with U\ having the op- 
posite sign than for the U\ used in the step simulations, 
because this sign of U\ gives the correct corrugation po- 



tential minima. For the step simulations it did not really 
matter what sign we chose, but for the point defect prob- 
lem it is more important to place the point defects at the 
correct positions. For example if we want Vi to represent 
the potential due to a vacancy we can choose Vq posi- 
tive and equal to the depth of the corrugation potential 
minimum, so that that minimum, presumed to be due 
to a substrate ion located at that position, is canceled 
out by Vq.) Because of the use of periodic boundary 
conditions, it was necessary to include the eight nearest 
neighbor cells to the cell containing the film in which the 
simulations were done. Each of these neighboring cells 
contained the defect potentials reflected into these cells. 
The impurity positions and the reflected impurity posi- 
tions are shown in Fig. |[ Only interactions between a 
film atom and a defect (or an image of a defect) which 
is less than 4er away from the atom are included in the 
calculation of the force or potential acting on the atom. 

Defect Positions 



(Angstroms) 



FIG. 3. The defect positions are shown. The central box 
is the cell in which the simulations were performed, and the 
neighboring cells represent the reflected defect positions. 



B. Results of the Simulations for the Step 

In most of our simulations we have used a corrugation 
strength (defined as 9 times the value of fU\{z) evalu- 
ated at the value of z at which Uq(z) is minimum) of 
2.025 meV ( f=l), the value that was used in Ref. [Q]. In 
order to examine the effect of using a larger corrugation, 
we have, in addition, done some runs at a corrugation of 
3.56 meV (f=1.69). 

We apply an external force to all the particles in the 
monolayer at the following angles with respect to the x 
axis (which is perpendicular to the step): 0° (force per- 
pendicular to the step), 16°, 37°, 58°, 66° and 90° (force 
along the step). The case in which the particles move on 
a substrate which is free of steps ( with the force along 
the x-axis) is also studied to make a comparison with the 
simulations that were done with the step present. 

In Fig. [|, the velocity is shown as a function of time for 
the various values of the applied force for the case of a no 
step present, and in Fig. || for the case of a step present. 
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The velocity at the top of the plateau is the steady-state 
velocity, which we will later plot versus the applied force 
F. In order to diminish thermal fluctuations of the ve- 
locity vs. time curves, the curves shown in Figs. ^ and 
^ are obtained by making averages of several runs. The 
no-step case was obtained averaging 5 runs and the step 
cases, by averaging 3 runs each. The steady-state veloc- 
ity is plotted as a function of F in Fig. [g. We see that 
the magnitudes of the velocities depend on the orienta- 
tion of F. The velocity is smallest at 0° (i.e. when the 
particles move almost perpendicular to the step) . It gets 
larger when the angle increases, taking on its maximum 
value for an angle of 90° (when particles move along the 
step). 



FIG. 5. Velocity of the center of mass in the x direction as 
a function of time, for various values of an applied force, for 
the case of a step present. In this case the driving force is 
applied at 0° with the x axis. 




Force (e/c) 



When the particles move perpendicular to the step 
(i.e., in the x direction), their motion in x is pinned. 
For angles of incidence < < 90 the pinning transition 
appears to be "softened" or "smeared out." It is thus 
of interest to plot the x and y components of velocity 
vs. the applied force (i.e. V£™ vs. F) in order to gain 
insight into what the atoms do. In Fig. 0, the x and y 
components of the velocity are shown separately for 58° 
for a 2.025 meV corrugation potential strength. We see 
evidence of pinning in x but not in y. (The exponent (3 
in this case is j3 = 1.4.) 
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FIG. 4. The velocity is shown as a function of time for 
the various values of the applied force for the case of no step 
present. The driving force is applied at 0° with the x axis. 
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FIG. 6. Variation of the velocity with F at 77.4° kelvin. 

[3=1.4 



Correlation coefficient = 0.9807333 

0.40 | , 1 , 1 , r- 




Force (e/o) 



FIG. 7. Variation of v x and v y with F for 58° for 2.025 
meV of corrugation. 

We also made some runs with a corrugation strength of 
3.56 meV. The film moves in the same fashion as for 
2.025meV corrugation, but of course with much smaller 
velocities. 

Two cases clearly show linearity between velocity and 
F: no step and 90° (along the step). In addition, for an- 
gles between 0° and 90° v shown in Fig. || is linear in F 
to a good approximation for v less than 0.1a/ to . This is 
one possible explanation for why such viscous friction is 
found in experiment. 

C. Stick Slip and Creep 

The results reported in Fig. || are obtained by averag- 
ing the V cm versus time data for a few different runs for 
each value of the applied force. These averages are per- 
formed in order to suppress the large fluctuations in the 
data found when a step is present. In order to illustrate 
this, V c m is plotted as a function of time, without doing 
any averaging, for one value of the applied force larger 
than the critical force in Fig. |[ For comparison, a plot 
of Vcm versus time for the case of no step present is also 
made on the same graph for an applied force chosen so 
that V cm has about the same value at each time. In this 
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graph the force, for the case in which there is a step, has 
a value F — 1.9 1CP 3 e/er, slightly above threshold. 



Across step F=1 .9 1 s/a 




time (s) 

FIG. 8. Comparison of velocity profiles for the step and 
no step case. Fluctuations in V cm are noticeably larger when 
there is a step than in the absence of it. 



It is clear from this figure that the the fluctuations in 
V cm are noticeably larger when there is a step than in 
the absence of a step. It is likely that these large fluc- 
tuations are due to stick-slip motion, which occurs for a 
system which would be pinned at smaller values of the 
applied critical force. 



In Fig. ^, a typical plot of V cm is shown for F normal 
to the step below its threshold value. The x component 
(normal to the step) and the y-component of V cm are 
shown separately. V cm>x becomes thermally activated for 
a short period of time and then becomes pinned again. 
Such behavior is not found in V cm . y . The behavior for 
values of the force slightly above threshold is shown in 
Fig. [h]. In Fig. [ll] a similar plot is shown for the applied 
force at an angle of 37 degrees with the normal to the 
step. We can again see the thermally activated behavior 
of V cm ,x for a couple of short time intervals. In contrast, 
V cm .y appears to saturate at a positive value (because 
the film is not pinned in the y-direction) . Thermally ac- 
tivated motion of the type that we see here is the type of 
behavior that would lead to creep of a macroscopic film. 
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FIG. 10. V C m vs. time for F > F c across the step. 
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FIG. 11. The applied force is at an angle of 37 degrees with 
the normal to the step. We can see the thermally activated 
behavior of V cm ,x for a couple of short time intervals. In con- 
trast, Vcm.y appears to saturate at a positive value (because 
the film is not pinned in the {/-direction) . Thermally activated 
motion of the type that we see here is the type of behavior 
that would lead to creep of a macroscopic film. 



D. Film Rotated Relative to the Step 



Across step F=1.5 1(T e/o 




FIG. 9. Typical plot of Vcm vs time for F < F c The x 
component (perpendicular to the step) and the {/-component 
of Vcm are shown separately. V cm ,x becomes thermally ac- 
tivated for a short period of time and then becomes pinned 
again. Such behavior is not found in V cm ,y 



In the runs discussed so far, the crystallographic axes 
of the film were lined up with the step. This is not al- 
ways true for this system. For example, we also did a 
run for a 173 atom film (a more compressed monolayer 
film) for the same substrate containing a step. In this 
case, the film in equilibrium APPEARED rotated at an 
angle with respect to the step, as seen in Fig. [jj]. In this 
case, the film attained a steady-state velocity of 0.2<r/to 
when F=0.0015e/fi was applied normal to the step. As 
can be seen in Fig. [g, the 163 atom film, which was not 
rotated with respect to the step, was pinned for F of that 
magnitude. As we shall see in the next section, this is 
consistent with the predictions of lowest order perturba- 
tion theory. 
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QCM friction measurements can be explained by rela- 
tively small value that we have found for the pinning 
force. We do not believe that this is correct, however, 
because this force is still much larger than the effective 
inertial force due to the oscillations of the QCM's sub- 
strate which is equal to muP'A m 10~ 14 <iyn, where m 
is a film atom mass (« 10~ 22 dyn, to « 10 7 s _1 and the 
amplitude of the substrate oscillations A w 10 2 A°. 



FIG. 12. The atomic positions (looking down on the sub- 
strate) are shown for 173 atoms (171 atoms adsorbed on the 
substrate and two remain above the substrate). A side view 
of the step is also shown at the bottom of the figure as a guide 
to the eye. This film is rotated slightly with respect to the 
walls of the box, and hence with respect to the step. 



E. Results of the Simulations done with Point 
Defects Present 
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In these simulations Vq in Eq. |2.6| was chosen to be 
equal to 0.1 e, for which each defect potential has a max- 
imum value approximately equal to the corrugation po- 
tential well depth, and the concentration c was chosen to 
be equal to 0.05. Again, the system was first equilibrated 
for at least 600,000 iterations, and then an external force 
was applied. Simulations were done with the external 
force making angles of 33, 105 and 203 degrees with re- 
spect to the side L x of the box. The results are shown 
in Fig. |l3|. Each point represents the average of the 
component of the center of mass velocity along the ap- 
plied force over a time interval equal to at least IOHq (i.e. 
2 x 10 6 iterations or about 3 x 10~ 8 s.) in order to aver- 
age out the fluctuations. An average of three runs is also 
shown for one of the cases. Since it did not look qualita- 
tively different than the single runs, it is clear that single 
runs are adequate. For all three directions of the applied 
force the film appeared to become pinned below a criti- 
cal force of about 2.10 -4 e/a, even though the detailed 
shape of the velocity versus force curves might have dif- 
fered slightly. A run was also done at a temperature of 26 
K in order to determine whether the motion of the film 
might have been partly thermally activated. The results 
are shown in Fig. |lj. It is seen that the critical field is 
now about 4.10 -4 e/a, and the curve appears to exhibit 
a sharp depinning transition, whereas the runs done at 
T = 77 A °K showed what looked more like a rounded 
transition. These results imply that some thermal acti- 
vation of the atoms out of the impurity potential wells 
is taking place. It is quite unlikely, however, that such 
thermal activation, which is responsible for creep, will 
give such a large contribution for a macroscopic film as 
it does for the small films used in the simulations. 

The pinning force of 2.10" 4 e/a found at T = 77.4°^ 
is equal to about 1.5 x 10 _11 dj/n per atom. 

One might be tempted to think that the fact that point 
defects do not seem to have much of an influence on 
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FIG. 13. The velocity in units of a/to versus applied force 
f in units of e/a is shown for a lattice containing a random ar- 
ray of point defects with a 5 percent concentration, the force 
f applied at 33, 105 and 203 degrees with respect to L x . 



An average of three runs is also shown for one case. 
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FIG. 14. The velocity in units of a /to versus applied force 
f in units of e/a is shown for a lattice containing a random 
array of point defects with a 5 percent concentration, with 
the f making a 33 degree angle with L x . This run was done 
at a temperature of T = 26 °K . 



III. PERTURBATION THEORETIC 
TREATMENT 

The force of friction acting on a thin film as it slides 
over a substrate containing defects was recently studied 
using lowest order perturbation theory in the substrate 
potential 0. It was shown there that for most orienta- 
tions of the film and most directions of sliding over the 
substrate, the force of friction is viscous (i.e., propor- 
tional to the sliding velocity), in agreement with what 
has to date been observed in QCM friction experiments 
||. The simulations that were discussed in the last sec- 
tion of a xenon film sliding over a silver (111) substrate 
containing a step show that the film is pinned if we apply 
a force smaller than a critical value in a direction normal 
to the step. When a force larger than this critical value 
is applied, the resulting velocity is a nonlinear function 
of the applied force. In this case the film's axes are lined 
up with the step. For many substrate-adsorbed film sys- 
tems, the case considered in the simulations, in which 
the crystallographic axes of the film are lined up with 
the step, is an important case because one expects that 
it will be energetically favorable for the film 's axes to line 
up with the step locally. The reason for expecting this 
to occur is that adsorbate atoms can lower their energies 
by lying along a step edge, because in this way each of 
these adsorbate atoms will be surrounded on two sides 
by a substrate atom. For the present system, our sim- 
ulations seem to show, however, that the film does not 
always line up with the step, and when it is not lined up 
(for example, for the 173 atom system discussed at the 
end of the last section), the simulations show that the 
film is not pinned for forces for which it would be pinned 
if it were. (For real surfaces, which have many steps, it is 
expected that because of its stiffness the film in its solid 
phase will typically not be able to distort in such a way 
that its crystallographic axes line up locally with each 
step.) In this section we will try to interpret the results 
of the simulations reported here using the perturbation 
theoretic methods of reference. j?]] 

In Ref . |?J , it was shown using perturbation theory that 
for a general nonperiodic substrate force, the mean force 
of friction acting on a thin film is given by 



F av v = (mN)- 1 J2 



7 o; 2 |f(k).eV| 2 
(w£(k) - uj 2 ) 2 + 7 2 w 2 



(3.1) 



where u> = v x k x + v y k y , where v is the sliding velocity of 
the film, f (k) is the Fourier transform of the force due to 
the substrate, m is the mass of an atom, N is the number 
of atoms in the film and 7 is the inverse phonon lifetime. 



This expression was obtained by setting the rate at which 
the substrate force f(f— vt) does work on the film equal 
to the rate at which the average force of friction F av does 
work on the film, i.e., F av v. It was also found in that ref- 
erence that the mean square vibrational displacement of 
an atom in the film (in the rest frame of the film) is given 
by 



< u >= m 



2 ^E 



|f(k) • e k , CT | 



(u> 2 (k) - LO 2 ) 2 + 7 2 W 2 



(3.2) 



In Ref. [ ?J , it was found that in two or fewer dimensions 
for a substrate containing point defects (e.g., vacancies 
or substitutional impurities), < u 2 > diverges as v ap- 
proaches zero. This signifies that in the zero velocity 
limit, there will always be significant distortion of the 
film to conform to the substrate, no matter how small 
the substrate force, which implies that the film will be 
pinned in place unless a strong enough external force is 
applied. 



In the 7 approaches zero limit, Eq. (3.1) becomes 



F av = (Tr/m)(a/27r) 2 J2 J d d ku\i (k)| 2 <5(^(k) 



(3.3) 



In this limit, a nonzero contribution to F av occurs when 
the argument of the delta function is zero, which occurs 
when the plane in the d+1 dimensional space defined by 
u> and the d-components of k whose equation is lu = v ■ k 
intersects the phonon dispersion surface [whose equation 
is uj = uj a (k)], if f(r) were aperiodic for all directions of 
r. This is the case when we have point defects in the 
substrate. It was shown in Ref. M that in this case in 
lowest order perturbation theory for a two dimensional 
(i.e., monolayer) film F av is independent of velocity and 
< u 2 > diverges as 1/v as v approaches zero, implying 
pinning of the film at low velocities. In contrast for a 
three dimensional (i.e., thick) sliding film, F av was found 
to be proportional to v and < u 2 > did not diverge as v 
became zero, which is consistent with the film not being 
pinned. In appendix C, these ideas are expanded upon 
by considering general order in perturbation theory. It 
is shown that d=2 is indeed a critical dimension for this 
problem. 

Let us now consider a line defect, such as a step or 
facet boundary. If we choose the x and y axes (i.e., the 
coordinate axes in the plane of the film) so that the y-axis 
is along the defect (i.e., along a step or facet boundary, 
which are taken in the present discussion to be straight 
lines) f (k) will be non-zero only if k y is a multiple of the 
y-component of one of the reciprocal lattice vectors of the 
substrate, since the substrate is periodic in that direction 
if the crystallographic axes of the substrate are lined up 
with the defect and quasiperiodic otherwise. As a conse- 
quence, the plane whose equation is uj — v x k x + v y k v gets 
fragmented into a series of parallel lines at values of k y 
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equal to Q v , the y-component of one of the smaller recip- 
rocal lattice vectors of the substrate. For a facet bound- 
ary or step, f can be written as f(r) = g(x)F(r), where 
F is periodic but g is not. (This was illustrated for the 
model for the potential that was used in the simulations 
presented in Ref. jjj). A reasonable model potential for 
other line defects, such as facet boundaries, for example, 
will have a similar form. This is discussed in appendix 
B. Since we can write F as F(r) = Fqe'Q r , we find 
that 



not pinned. Furthermore, if |f • el is of the order of 



f(k) = / 



d 2 , 



Q 



,Q y ffk„-Q.FQ. ( 3 - 4 ) 



where Q denotes a reciprocal lattice vector of the sub- 
strate. It is only when one of these lines in k-s pacc inter- 
sects the phonon dispersion surface that Eq. (3.3) gives 
a nonzero contribution to F av . If this does not occur, we 



must keep 7 nonzero in Eq. (3.1). We then find that F av 
is proportional to v |Q . 

Let us now consider the case in which the crystallo- 
graphic axes of the film are not lined up with the direction 
of the defect but are close to being so. Furthermore, let 
us keep 7 nonzero for this discussion, in order to examine 
what happens if a line passes within a phonon linewidth 
of the phonon dispersion surface Then taking the defect 
to be lined up with the y-axis . if we convert the summa- 
tion over wavevector in Eq. (3.2) to an integral along a 
line perpendicular to the step in the usual way, we obtain 



< u 2 > /a 2 



2nm 2 a 2 ^ 



dk x F(k x ,k y ,G,v) (3.5) 



where F(k x , k y , G, v) is given by 
\{{k x , k y0 ) ■ e k , CT | 2 



[^(^ + ^o)-^ 2 ] 2 +7 2 w 2 ' 



(3.6) 



where to — v y k y Q + G x v x and where k y Q is the amount 
that the line, which would have passed directly through 
the point in k-space denoted by the film reciprocal lattice 
vector G if the film's crystallographic axes were line up 
with the defect, misses going through this point in recip- 
rocal space. On doing the integral over k x by contour 
integration we find that < u 2 > /a 2 is equal to 



/ \ — 2 — 1 / \ 
(ma) v (— ) 



7r|f • 



2^'2 Sl n(e 1 /2)(v^ y0 + 1 2 G 2 x vlfl^ 



(3.7) 



where |f • e| is evaluated at the wavevector correspond- 
ing to the pole in the integrand of Eq. ( |3.5| ) and where 
9 1 = arctan[yv x G x / \v 2 G 2 x — v 2 k y0 )]. (Note that sin6i/2 
becomes 1 rather than zero as v x approaches zero, since 
9\ approaches 7r rather than zero. ^From these results, 
we see that if k y o 7^ 0, < u 2 > /a 2 does not diverge 
as v x approaches zero, which implies that the film is 



10- 7 dyn, 



m 



10 



-22, 



10" s cm, G v 



10 8 cm" 



and k y o ss 10 cm (or about 1/10 of the Brillouin zone 
radius), we estimate that < u 2 > /a 2 is of the order of 
10~ 3 . Thus, we conclude that for a reasonable value of the 
mis-orientation parameter k y o , lowest order perturbation 
theory should be a correct description of the friction, even 
as v approaches zero. It is easily seen from Eqs. (3.1) 
and (3.2) that for this case, in which < u 2 > /a 2 does 
not diverge as the sliding velocity approaches zero, the 
force of friction will be proportional to the velocity. 



IV. SPECULATIONS ON THE BEHAVIOR OF 
REAL INTERFACES 



We propose the following physical explanation for why 
a line defect, such as a step, will pin the film if it is ori- 
ented along a crystallographic axis of the film but will not 
pin it if it is not lined up: When the line defect is lined 
up with a crystal axis, a line of atoms running across the 
width of the film will pass over the defect all at the same 
time as the film slides over the defect. Consequently, the 
force exerted by the defect on the film will be an extensive 
quantity, in the sense that it will scale with the width of 
the film. In contrast, when the defect is not lined up with 
an axis, atoms in the film will pass over the defect one at 
a time. The defect will act on one atom at a time with 
a force that does not scale with the width of the film. 
Thus, if we consider a situation in which the velocity of 
the film is studied as a function of a force applied to each 
atom in the film, when the film's axes and the defect are 
not lined up, it will clearly be much easier for the applied 
force to overcome any nonextensive pinning force due to 
the defect than an extensive one. A real substrate sur- 
face will contain a finite density of line defects along a 
line drawn along the sliding velocity. These defects will 
most likely not run across the width of the film, and they 
certainly will not be straight over that distance. Then, 
although the pinning force due to one of these line defects 
will not be extensive, the force due to a finite density of 
them will. 

One way to understand the case of many finite length 
line defects, which are not perfectly straight, as occur on 
real surfaces using the calculation presented above for a 
single straight line defect is that if the defect is not a 
straight line over an infinitely long distance , th e Fourier 
transform of the substrate force in Eq. (3.4) will no 
longer contain a 6k ,q factor (since a kroniker delta only 
occurs if f is periodic over an infinite distance). Rather 
the Kroniker delta will be replaced by a function of k y of 
finite width of the order of 2ir/£, where I is the length of 
a straight section of the defect, peaked around the values 



of Q y . Consequently, we can now think of the lines in 
k-spacc, defined by k y — Q y , which must intersect the 
phonon dis persion surface in order for the delta function 
in Eq. (3.3) to be satisfied, as being broadened. For real 
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surfaces, which contain many hnite length straight sec- 
tions of line defects separated by defect- free regions, f (k) 
will consist of the following contributions: defect-free re- 
gions which will give a peak in f at each reciprocal lattice 
vector of the film, point defects, which will give a contri- 
bution which is not peaked in k-space and straight sec- 
tions of line defects, which will each give a contribution 
that is peaked at values of the component of k along the 
defect equal to components of a reciprocal lattice vector 
of the substrate if we move in k-space in a direction par- 
allel to the defect. How sharply it is peaked will depend 
on the length of the straight section of the defect. In the 
direction perpendicular to the defect, the contribution to 
f (k) will not be peaked. The magnitude of each of these 
contributions to f (k) will be proportional to the number 
of lattice sites in each of the above elements. Th us, when 
we t ake the square of f(k), which enters Eqs. ( |3.l|) and 
(3.2), we will thus obtain the following contributions: 1. 
a Bragg peak at each substrate reciprocal lattice vector, 
2. a contribution which is not peaked in k (from the point 
defects), and 3. a series of "mountain ridges" in k-space 
directed in various directions. These are the lines in k- 
space that we found for a single line defect except that 
each one now has a width of the order of 2ix divided by 
the length of the straight section of the line defect that 
gives this contribution. If the widths of the "mountain 
ridges" are fairly narrow, the ridges will generally miss 
passing through a dip in the phonon dispersion surface, 
which as discussed in the la st se ction, is the criterion for 
the delta function in Eq. (|3.3|) being satisfied). Then, 
these will give mainly a viscous contribution to the fric- 
tion, following the discussion in the last section. On the 
other hand, if the straight sections of the line defects 
are relatively short, resulting in broad ridges, the ridges 
will generally intersect the phonon dispersion surfaces, 
most likely leading to "dry" (i.e., velocity independent) 
friction. The point defects will always contribute "dry 
friction." We know from the discussion in the last section 
that the Bragg peaks will give viscous friction in pertur- 
bation theory, as long as their widths are not too large. 
More numerical studies of this kind, which use pertur- 
bation theory to study the effects of realistic models for 
defected surfaces will be presented in future work. 



V. CONCLUSIONS 



Simulations and perturbation theory calculations done 
for a monolayer film sliding on a substrate containing 
a step defect gave approximately viscous friction for all 
sliding directions other than perpendicular to the step. 
The point defects in the simulations were found to pin the 
film below a critical applied force For defects of strength 
comparable to the strength of the corrugation potential 
(which were used in the present simulations), however, 
the pinning force was only about 10~ 11 dyn. This force 
is still greater than the inertial effective force resulting 



from the oscillations of the microbalance (mw 2 4, where 
m w 10 _22 <7 is the adsorbate atom mass and to « 10 7 s _1 
and A w 10 _6 cm are the microbalance frequency and 
amplitude, respectively), which is about lO^^dyn. The 
fact that the observed force of friction generally seems to 
be viscous in these experiments seems to imply that ei- 
ther the defects are much weaker or have a concentration 
much lower than the 5 percent concentration used in the 
present simulations. 
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APPENDIX A: EFFECTS OF PERIODIC 
BOUNDARY CONDITIONS 

We will now use perturbation theory Q to estimate 
the effect of using periodic boundary conditions in sim- 
ulations of films sliding over substrates containing steps 
and substrates containing point defects (e.g., vacancies, 
substitutional impurities or add-atoms). What we will 
do is to calculate the force of friction for a system with 
periodic boundary conditions and compare our results to 
the results of the calculation of the force of friction for 
an infinite system. For the case that we consider, mo- 
tion takes place perpendicular to the step. Since the fi- 
nite system with periodic boundary conditions is periodic 
with the length of the box, the driving force due to the 
substrate (in which we are doing perturbation theory) is 
periodic. As a result, the sum over wave- vector in the ex- 
pression for the force of friction, i.e., Eq. (5) in Ref. 0, 
becomes Eq. (12) in that reference, which gives the force 
of friction when the substrate is perfectly periodic with 
reciprocal lattice vectors with components G x = 2im\ / L 
and G y = 2irri2/L, where L is the length of the box (the 
reciprocal lattice vectors for a substrate which is peri- 
odic only because of the use of periodic boundary condi- 
tions). Then the important parameters are R, the ratio 
of 2w/L to the reciprocal lattice vector of the substrate, 
the phonon mode width 7, and the sliding velocity v. 

To be explicit, we do the sum over wavevector in Eq. 
(5) of Ref. 0. When periodic boundary conditions are 
used with a finite size box, this is a sum over discrete val- 
ues of k (i.e., k x = n\2-K I L and k y = ni2-K / L, where L 
is the box length). We make the following assumptions, 
valid for small sliding velocity: We take the sliding ve- 
locity to be along the x-axis. We assume that the major 
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contribution to this summation comes from the intersec- 
tion of the plane uj — vk x and uj — ui a (k) . In the small v 
limit for point defects, Eq. (5) of Ref. can be approx- 
imated by 



F a v OC ^ 



(v 2 k 2 - v 2 G 2 x ) 2 + 7 2 w 2 



(Al) 



where G x is the x-component of one of the smallest re- 
ciprocal lattice vectors of the film. The k-dependcncc 
of f(k) • £k,o- is neglected and the origin of k-space has 
been moved to the point G, one of the smallest recipro- 
cal lattice vectors of the film. Then, we must evaluate 
the summation 



E 



rv 



[{n\ + n 2 )R 2 



,'212 



(IV) 



(A2) 



where R = 2tt/(LG x ), v' — v/v p and T ~ ^y/vpGx. This 
summation is evaluated for R=l/13, T = 0.1 and v' rang- 
ing from to 0.2, which are parameters appropriate to 
the simulations. In the infinite L limit, the summation 
over ni and ri2 can be approximated by the integral 



where q x = 
integral is 



dq x dq y 



Tv' 



(q 2 - v 12 ) 2 + (TV 1 ) 



A2 ■ 



(A3) 



n\R and qi — n^R. The result of doing this 



(R- 1 )[n/2 + arctan(v' /T)]. 



(A4) 



Plots of Eqs. flA2| ) and ( |A4j ) are shown in Fig. |15|. As can 
be seen, the integral (which describes an infinite system) 
and the summation (which describes a finite system with 
periodic boundary conditions) agree with each other for 
v' greater than 0.02. Thus, for velocities greater than 
this value, we are justified in assuming that the simula- 
tions will be a good description of a film on an infinite 
nonperiodic substrate. 




0.00 0.05 0.10 0.15 0.20 

v' 



FIG. 15. A comparis on o f Eq. (A2) (the curve that is 
higher near v'=0) and (A4), which are proportional to the 
force of friction for a finite and an infinite system, respec- 
tively, with periodic boundary conditions. 



For the case of a straight step, the film is periodic or 
quasiperiodic along the step with a repeat distance or 
almost repeat distance comparable to a lattice constant. 
Thus, if we take the step to be along the y-axis and con- 
sider the sliding velocity to be in the x-direction, the 
expression for F av involves an integral only over the x- 
component of the wavevector for an infinite system. For 
a finite system with periodic boundary conditions, the 
summation and integral corresponding to these quanti- 
ties are respectively 



E 



Tv' 



(n 2 R 2 - v 12 ) 2 + [Tv'f 



and 



R~ 



dq- 



Tv' 



(A5) 



(A6) 



(q 2 - v' 2 ) 2 + (Tv') 2 ' 

where q=nR. This integral can easily be performed by 
contour integration to give (n/ (27) 1 / 2 (u')~ 1 / 2 . This re- 
sult together with the summation over n are shown in 
Fig. [U]. We again see that the finite system with pe- 
riodic boundary conditions is an accurate representation 
of the infinite system for v' greater than 0.02. 
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V' 



FIG. 16. A comparis on o f Eq. (A5) (the curve that is 
higher near v'=0) and (A6), which are proportional to the 
force of friction for a finite and an infinite system, respec- 
tively, with periodic boundary conditions. 



APPENDIX B: THE SUBSTRATE FORCE WITH 
A LINE DEFECT PRESENT 

A good model for the corrugation potential when there 
is a step present in the crystal which is parallel to the y- 
axis is 



vi(z - 



<?(*)) E' 

Qi 



,iQ 1 -(r-d(x)) 



(Bl) 



where g(x) is a Fermi function z\l\e^ iiyX ~^ + 1] times z\ 
(where Z\ is the height of the step and (3 is the reciprocal 
width of the step edge, which is of the order of the recip- 
rocal of an atomic spacing), Qi denotes the 6 smallest 
reciprocal lattice vectors of the substrate and d(a;) is the 
product of a vector giving projection of the difference be- 
tween the equilibrium positions of the film atoms in the 
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first and second atomic layers parallel to the x-y plane 
and the above Fermi function. Then the x and y compo- 
nents of the force f x / y due to the corrugation potential 
are given by 



Z Qi X 



Qi 



iQi-(r-d(x)) 



and 



fy = Vx 



Qi 



(r-d(x)) 



(B2) 



(B3) 



where zq is the value of z at which the attractive poten- 
tial Vq(z) is minimum. Then it is easily foun d th at each 
term in f x and f y will have the form of Eq. (3.4). 

Now, let us consider another type of line defect, a facet 
boundary. In going from one facet to its neighbor, the 
crystal surface rotates in the x-z plane for a facet bound- 
ary running along the y-axis. Then a simple phenomcno- 
logical model for a facet boundary, along the lines of the 
model discussed above for the step, is obtained by allow- 
ing the rotation angle to depend on x. Then, we may 
write 



x' = xcos6(x) — zsin9(x) w x — zO(x), 



(B4) 



z = zcosQ(x) + xsin8(x) z + x9(x), (B5) 
since the rotation angle 9(x) is small. Then we obtain 



fx = Y}- dV ^ dz '\z'=z Q {0{x) + x6' (x)) + 
Qi 

»Qix^i(«b)(l- z6'{x))]e 



?i«*of(aO-QiT 



whose Fourier transform will clearly be proportional to 
_> ly , as is Eq. (3.4). The fundamental difference be- 



tween these two line defects is that whereas dg/dx, which 
determines the strength of the nonperiodic part of the 
substrate force (due to the presence of the step) is of 
order 1, 9, which determines the strength of the nonperi- 
odic part of the substrate force due to the facet boundary, 
is much smaller, since 6 is the angle between two neigh- 
boring facets, which is likely to be only about a tenth of 
a radian. Thus the force due to a facet boundary is an 
order of magnitude smaller than that due to a step. 



APPENDIX C: BEHAVIOR OF A GENERAL 
TERM IN PERTURBATION THEORY 



examine the behavior of the n th order term in perturba- 
tion theory in the substrate potential as the sliding veloc- 
ity approaches zero. (The charge density wave problem 
differs from the present problem, however, in that in the 
charge density wave problem, the motion is overdamped, 
whereas in the present problem it is not.) To accomplish 
this, let us consider the average rate at which the sub- 
strate force does work on the film by generating phonons, 
as the film slides over it (which can be set equal to F av v 
where F av is the mean force of friction and v is the slid- 
ing velocity). From previous work on this problem this 
is given by |J 

F av v = (T- 1 ) I dt f (Rj + Uj - vt) ■ Uj(t), (CI) 
J j 

where 

«i(*) = E / dt'GiRj-R^t-t')- 

j' 

f(R r +n jl (t')-vt'), (C2) 
where the Green's function is given by 

EC e tk-(%-Ity) e u>(t-t') 
/ (hj^—r ; . 

(C3) 

The perturbation series in the substrate potential corru- 
gation strength f can be generated by ex pan ding f (Rj' + 
Uj' — vt') on the right hand side of Eq. ( |C2j ) in a Taylor 
series in Uj/(i'). Each factor of Ujr in each term in the 
Taylor expansion is given by Eq. (C2). We then make 



a Taylor series expansion of the f (R jV + Uj'(i') — vt') 
which occurs in each of these expressions. At any given 

point in this procedure, we can stop it by replacing 
(B6) f(Rj,+Uj/(t')-v£')byits 

zcroth order term in the Taylor 
series expansion, namely f (Rj' — vt') . Then the resulting 
perturbation expansion of Eq. (C2) can be represented 
schematically by the diagram in Fig. [l?]. 




Arguments similar to those used by Fisher for a similar 
problem, that of a charge density wave moving in a crys- 
tal lattice in an applied electric field |10|, can be used to 



FIG. 17. Diagramatic representation of the perturbation 
expansion for Uj in powers of f. The lines and vertices are 
defined in the text. 
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Each vertex is labeled by a position and time (Rj n ,t n ), 
where n is equal to 1 at the vertex that is furthest to the 
left. It increases as one moves to the right. With the ex- 
ception of the vertex at the beginning of the "tree" (i.e., 
the far left hand side of the diagram) each of these vari- 
ables is summed and integrated over. A line connecting 
two vertices labeled by (Rj 1 ,t\) and (Rj 2 ,t2), for exam- 
ple, represents the Green's function G(R J1 — Rj 2 , t\ — tz), 
and each vertex denotes (q — 2) nd order derivative of f 
multiplied by a factor of l/(q — 2)!, where q is the num- 
ber of lines intersecting at the vertex. It is not difficult 
to show that the last element in a chain of vertices (i.e., 
a line of vertices, each one of which has only two neigh- 
bors), which ends at the vertex (Rj^ii) is a function 
of Ri and t\ only in the combination (Rj 1 — vii). For 
example, in d-dimensions this element is proportional to 



that a chain ending at the vertex labeled by (Rj i; ti) is 
a function of the variable (R^ — vfi), it is easy to show 
that a vertex at which q chains intersect will also have 
this dependence on the space and time variables with 
which it is labeled, and it will also contain a factor 



IEU> 2 (k m ) - (v • k m ) 2 + i 7 k m )]-\ 



(C9) 



/ 



d d k 



;*-(Ri-vt) f ( k ) . £k 

Lu 2 (k) — (v • k) 2 + ijv ■ k 



(C4) 



In two and three dimensions the above integral can be 
done by contour integration. In the small v limit w(k) 
can be approximated by v p k', where v p is the sound ve- 
locity and k = G + k', where G is the nearest reciprocal 
lattice vector to the region in which w(k) and v • k inter- 
sect. Then the above integral becomes 



where n is the order in perturbation theory. Each variable 
k m is summed over. When we take the thermodynamic 
limit for the film, each of these summations becomes a 
d-dimensional integral over the wavevector. If these in- 
tegrals diverge at small v, the most divergent part will 
come from the integral over the vicinity of the intersec- 
tion of the plane w = v • k m with the phonon dispersion 
surface Lu(k m ). Therefore, let us consider the integral 
over just this region of k-space. Then in the small v 
limit, uj(k m ) can be expanded in a Taylor series about 
the film reciprocal lattice vector G of the dip in cu(k m ) 
which is intersected by the plane u> = v • k TO . Then, if we 
write k TO = G + k' m , v • k m and oj 2 (k m ) are, to lowest 
order in k' m , v • G and (v p k' m ) 2 , respectively. Then the 
divergent part of one of these integrals is proportional to 



/ 



d d k' 



v 2 k' [ 



(v • G) 2 + 17V • G ' 



(CIO) 



d d k' 



iMRi-vt)f( G ) . £G 



(C5) 



Vpk' 2 - (v • G) 2 + 27V • G ' 
In one dimension this integral is proportional to 

cos(vG(x 1 - vh))e-^ Xl - vtl \. (C6) 



Let us consider this integral in the limit as 7 approaches 
zero since it is expected that the phonon damping will al- 
ways be small compared to the frequency of the phonon 
being excited (which is equal to v • G). If we make a 
change of integration variable to a variable x equal to 
v p k/(vG) the integral becomes 



In three dimensions, it is proportional to 

e »(vG/t;p)|Ri-vt| e -(7/t;p)|Ri-vti| 

i 1 • 

Ri - vii 



{vG/v p ) d {v p /vGf J 



x.(R 3 -vt) 

+ ie ' 



(Cll) 



(C7) 



By building a chain out of these elements and integrating 
over each of the vertices except the furthest to the left, 
we find that if (Rj^ti) are the labels of the end of the 
chain farthest to the left, the chain will be a function of 

these variable in the combination (R^ - vh). On do- finite number. Then the n th order term in the perturba 



where e = j/vG, which we are assuming to be << 1, 

since the damping of an acoustic phonon mode is always 

small compared to its frequency. The smallest value of 

x in the range of integration is and the largest value is 

large compared to 1, but still finite. As long as e is small 

compared to 1 but nonzero, the integral over x gives a 

th 



ing the integrals and sums on each of the vertices in the 
chain except the last one, the Fourier transform of the 
n th Green's function in the resulting contribution to the 
perturbation theory is 



tion theory diverges as v d ~ 2 if d < 2. Two dimensions 
is a marginal dimensionality, which could still have loga- 
rithmic divergences in v at small v. 



[ W 2 (k„)-(vk„) 2 +J7(vk„)]- 



(C8) 



where k„ is one of the wavevectors which are summed 
over. Each vertex represents the dot product of one of the 
phonon polarization vectors ek„ and the spatial Fourier 
transform of one of the functions f, evaluated at the dif- 
ference between the wavevectors of the Green's functions 
on either side of the vertex, multiplied by the dot product 
of a polarization vector and the argument of the Fourier 
transforms of f appearing in the chain. Using the fact 
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